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Abstract 

Hierarchy structure of the Bethe-Salpeter (BS) equations is found. We also 
employ this structure to greatly improve on the approximations for the BS 
(-H , kernel. Resummation of the BS kernel in t and u channels to infinite order 

CM is equivalent to truncate the effective action to infinite order. Our approxi- 

p [ mation approaches insure that the theory can be renormalized, which is very 

important for numerical calculations. Two-point function can also obtained 
from the four-point one through flow evolution equations resulting from the 
exact renormalization group. BS equations of different hierarchies and the 
flow evolution equation for the propagator constitute a closed system, which 
can be solved completely. 

Keywords: 

Bethe-Salpeter equations, Renormalization, Exact renormalization group 



1. Introduction 

Bethe-Salpeter (BS) equation plays an important role in many-body the- 
ories. It has been widely used in various research fields, from strongly corre- 
lated electron systems |2[ , hadron physics 

mm 

, to particle physics and 



3n syste: 

field theory 0, 11, @, llDl , and so on. The BS equation resums its kernel, usually 
called as the BS kernel, to infinite order. Therefore, non-perturbative effects 
are included in the BS equation. However, usually in actual applications of 
the BS equation, it is impossible to obtain an exact kernel and we have to 
approximate it. For example, when the BS equation is employed to study 
the properties of mesons in the QCD, the usually adopted approximation is 
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called the rainbow-ladder approximation, in which bare quark-gluon vertices 
are used to construct the four-quark kernel However, the rainbow-ladder 
approximation is the simplest approximation, which loses lots of information. 
Lots of efforts have been done to go beyond the rainbow-ladder approxima- 
tion 

In fact, when we employ the BS equation to describe a many-body system, 
the central point is whether the kernel used is good enough to obtain what 
are observed in experiments. If it is not enough, how can we improve it? 
In this work, we will discuss how to improve the approximations for the 
BS kernel systematically. It will be found that we can greatly improve the 
approximations by employing the hierarchy structure of the BS equation. 

In this paper we will work under the formalism of the two-particle ir- 
reducible (2PI) effective action theory which is also known as the 



derivable approximations in the condensed matter physics [12|, ll3[. Recent 
years, the 2PI effective action theory has attracted lots of attentions. It has 
been applied to calculate the entropy of the quark-gluon plasma and other 
thermodynamic quantities 



14, 15 



with subsequent late-time thermalization|l6 



describe the non-equilibrium dynamics 
compute the shear viscosity in 
the thermal field theory 10(|, and so on. Specially, we would like to emphasize 
that due to many people's contributions 0, 0, M, [13]) has been clear that 
the 2PI effective action theory can be renormalized, which is quite non-trivial 
for a non-perturbative approach. Therefore, in our following discussions, the 
renormalizability of our approximation approaches will work as a prerequisite 
demand. The renormalizability should not be violated at any case. In fact, 
it is the demand of the renormalizability that leads us to find the hierarchy 
structure of the BS equation. 

The paper is organized as follows. In section [2] we discuss the hierarchy 
structure of the BS equation and how to improve on the approximations 
for the BS kernel by employing this hierarchy structure. In section [3] we 
show how the propagator and the effective action are obtained from the 
four-point vertex throu gh flow evolution equations resulting from the exact 
renormalization group 
outlooks. 



ISl . Il9| . Sec. m summarizes the results and gives 
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2. Hierarchy Structure of the BS Equation 



We consider the following scalar field theory with a non-local regulator 
term 

S4ip] = S[ip] + AS4ip], (1) 

with 

= - (2) 

AS^[ip] = -^(piR^^ijipj, (3) 

where iGf^lj = (—9^ — rn?)5ij, Here summations or integrals are assumed for 
repeated indices. The non-local term A5'k in Eq. (IHl) is employed to suppress 
quantum fluctuations whose wave lengths are larger than some value (here is 
denoted as I/k and k has a dimension of momentum) [l8|. This is achieved 
as follows: the non-local term in the momentum space is given by 

= -\j £^RMv{q)^{-q)- (4) 

The regulator Rniq) is chosen to have following properties: when g ^ k, 
Rk{ci) ~ K,"^, then the non-local term becomes a mass term with large mass 
K, which suppresses quantum fluctuations with wave lengths 1/q 3> I/k; 
when q ^ 1^, Rk{<i) 0, so fluctuations with wave lengths 1/q ^ l/n are 
not affected. In the exact renormalization group theory, beginning from a 
classical action at an ultraviolet scale A, one can obtain the corresponding 
quantum action which takes into account all quantum fluctuations through 
the evolution of flow equations from k = A to k = 0. 

From Eq. ([1]) we can obtain the corresponding 2PI effective action. The 



generating functional with one- and two-point sources is given by 

Z^[J, J2] = J [df] exp |i(^S'«[v9] + JiLfi + ^v9i J2,ij</9j j |, (5) 

W^[J,J2] = -tin Z^[J,J2]. (6) 
Performing functional derivative with respect to sources, one obtains 

-jj- = {Vi) = (7) 

1^ = ^(M + G.,), (8) 
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where is the expected value of the field Lp and G is the full propagator. The 
2PI effective action is obtained through the Legendre transformation from 
W^[J,J2l i.e., 

r^[0, G] = W^- Jt^y- - J2,ij- 



SJi ' 5J2,ij 
- Ji<pi - ^J2,ij{(pz4>j + Gij). (9) 



Expressed in terms of and G, the 2PI effective action is given by [11| 



r.[0,G] = ^<PiiGol,,<Pj+'-TTlnG-'+'-TTG,^G + r,^,[<p,G], (10) 



where we have ^Gqk ~ ^ ~ and the interacting part of the effective 
action is 



X + 5X 2 



2 u 4, ^ 

G + r2[0,G], (11) 

where r2 is given by all two-particle irreducible vacuum graphs whose vertices 
are given by the terms cubic or quartic in ip in the expanding expression 
of S[(f) + ip] — S[(j)], and propagators are the full ones. We have employed 
renormalized quantities in Eq. fIlOp and Eq. (ITT]) . They are related to the 
bare quantities (with subscript B) through the following relations: 

5m^ = Zml - , 6X = Z^Xb - A , 5Z = Z -1. (12) 
ZGq^ = Go 1 + 6Gq^ , 6Gq^ = i{6Zd^ + 6m^) , Gb = ZG. 

In actual calculations, we must make approximations. In another words, 
we have to truncate the two-particle irreducible vacuum graphs included in 
r2 in Eq. (ITT!) . For example, one can expand r2 in order of loop of skeletons, 
which is also known as the $ derivable approximation; one can also expand 
r2 in order of 1/A^ in the 0{N) model, and so on. Whatever it is, the 
common feature of these approximation approaches is that approximations 
are made to the effective action. Once the approximations are made, we 
can employ the approximative effective action to obtain the self-consistent 
equation for the two-point function. We can also obtain the BS equation for 
the four-point function. 
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In this paper, we will not go the same way as described above. We will re- 
verse the procedure, i.e., first, we make an approximation to the BS equation, 
then return to the two-point function and the effective potential. Comparing 
our approach with the conventional one, we find that our approach has several 
advantages: First, our approach makes much more powerful approximations. 
Second, we can easily observe the hierarchy structure of the BS equations 
from our approach. Finally, one can employ the hierarchy structure of the BS 
equations to reorganize the infinite diagrams through resummation. Same as 
the conventional approach, our approach produces approximations which are 
always renormalizable. As for the non-perturbative approximations, renor- 
malization is a quite non-trivial demand. 

In order to simplify our calculations but without lose of the generality, 
we consider the symmetric case in the following discussions, i.e., = and 
3-point vertex is also vanishing. We begin with the BS equation for the 
four-point vertex in the coordinate space: 

Mij-kl = ^ij;kl + -^-^ijik'l'Gk'i'Gl'j'Mi'j'.ki, (13) 

where the kernel is given by 

A.;M = 4.7^^%^- (14) 

We know that the BS equation resums diagrams in one channel to infinite 
order. Here we designate this channel as s channel. It is also known that the 
renormalizability of the BS equation demands that contributions to the kernel 
A only come from the t and u channels, and the s channel is prohibited in the 
kernel , [oj . As we have said, the s channel is resummed to infinite order 
through the BS equation. Therefore, it is natural to remind us a question: 
can we resum the two other channels to infinite order as well? In fact, this 
can be obtained by employing the BS equation once more, i.e., the kernel 
of the BS equation in Eq. f lT^ can be the solution of another BS equation. 
Therefore, we can express the kernel as 

A^r,kl = -^(A + 6X)6,,6,k6u + (M^,^^, - A^,^^,) + {M^j, - K^^^,), (15) 

and we have 

^il;jk = Ail;jk + 7^ AjZ;nj2 ^nii ^«2 j2 ^iii2;ifc- (1''') 
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It is easily checked that the kernel A has the same symmetry as A': Kij-^i = 

^ij;lk = ^ji;kl = ^kl;ij- 

It is more convenient to discuss the renormalization in the momentum 
space. In the momentum space, we attach momenta p, —p, —k, k with 
external legs whose subscripts are i, j, k, I, respectively. Then, the BS 
equations in Eq. flTB]) can be written as 

1 f d^q 

M'{p,-k;-p,k) = A'{p,-k]-p,k) + - J j^-:^A'{p,-k;-q,q + k-p) 

xG{q)G{q + k - p)M'{q,p - q - k; -p,k). (18) 

Here the kernel A' can be obtained by the corresponding 2PI vacuum di- 
agrams, same as that of the conventional BS equation, which also need to 
be truncated. For example, we can use 2PI vacuum diagrams up to two- 
loop or three-loop to obtain A', which corresponds to truncate A' to bare 
vertex or one-loop four-point function, respectively. Due to many people's 
efforts and contributions, we have known that the BS equation in Eq. f lTSj) 
can be renormalized, only when the kernel A' is two-particle irreducible, 
which isguaranteed by the two-particle irreducible property of the vacuum 
graphs U [s], lof. In order to make M' being finite, we need a count term in 
the kernel A'. Therefore, we can express A' as a sum of a finite part and a 
divergent constant, i.e., 

A'{p, -k; -p, k) = A'/p, -k; -p, k) - i6Xt, (19) 

where the subscript t indicate the t channel. In the same way, the BS equation 
in Eq. f|T7|) in the momentum space is 

1 f d^q 

M'{p,k;-p,-k) = A'{p,k]-p,-k) + - J j^-^A' {p,k]-q,q - k - p) 

xG{q)G{k + p- q)M'{q, k + p-q;-p, -k). (20) 

Also, we write the kernel as 

A'(p, k- -p, -k) = A'f{p, k; -p, -k) - idXu, (21) 

Employing 6X = 6Xs + SXt + 5A„, we express Eq. ( fTSl) in the momentum 
space as 

A{p,-p;-k,k) = -i{X + 6Xs) + M'{p,-k;-p,k) - A'j{p,-k;-p,k) 

+M'{p, k- -p, -k) - A'/p, k- -p, -k). (22) 
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One can see that A is also a sum of a finite part and a divergent constant 
(5As), wliicli make the four-point vertex M in the following BS equation 
finite: 



1 f d^q 

M{p, -p; -k, k) = A{p, -p; -k, + 2 J (2^^*^^' '^^ 

xG\q)M{q,-q;-k,k). (23) 

We should emphasize that diagrams included in the kernel A through resum- 
mation by Eq. ( JT8l) and Eq. (l20l) are two-particle irreducible in the s channel, 
which guarantees that Eq. fl23|) can be renormalized. 

Here we give some examples. The simplest case is that the kernel A' in 
Eq. ( fTSl) is just the bare vertex. Then M' is only dependent on p — k and 
Eq. ( |T8l) can be simplified to 



M'{p-k) = -t{X + 6Xt) + ^[-t{X + 6Xt)]M'{p-k) j 



xGiq)G{q + k-p). (24) 

Dividing both sides of the equation with — ^(A + 6Xt)M'{p — k) and using the 
free propagator Go, we obtain 

1 1 



-i(A + 5At) M'{p-k) 



' [\x\n— — ), (25) 



2(47r)^Ve Jq — x{1 — x){p — ky 

where we have employed the dimensional regularization and /i^ = 47r/i^e~'^-^ 
is a mass scale. If we introduce the following renormalization condition, 

' ' ' (26) 



-i{X + 6Xt) -iX 2(47r)2e 
i.e.. 



5X 



l-TjTT^- ^V2(47r)2e 

2(47r)2 t n=l ^ ' 

Then we have 



Finally, we get the kernel A as given by 
A(p, -p; -k, k) = i\ + 



^ ^ 2(47r)2 Jo "-^ m2-x(l-x)(p-fc)2 

+ ^ , ,1,7^ ,2 - ^SXs, (29) 

-'- + 2(47r)2 Jo (^^^^ m2-x(l-x)(p+fc)2 

One can easily find that the kernel A has the asymptotic behavior: A(p, —p; —k, 
A{p, —p; —k, fc) ~ at large k, which is the key point to renormalize the 
BS equation in Eq. ( l23l) . 



A' 





Figure 1: Contributions to the kernel A^j,.^-; up to one- loop. 



Then we consider a more complicated case and include one-loop contri- 
butions to the kernel A' in Eq. f|T6l) . which is shown in Fig. [H From the BS 
equation in Eq. ([T6|) . one can rewrite M' as 



M' = A' + ^A'G^A' + ^A'G^A'G^A' + (30) 

where the subscripts are not labeled explicitly. One can see that the second 
and third terms on the right hand side of the equation above correspond to 
one and two iterations. Substituting the kernel A' in Fig. [H we get their 
diagram representations which are shown in Fig. |2] and Fig. [3l respectively. 
One can observe that all diagrams in Fig. [2] and Fig. [3] are two-particle 
irreducible in the s channel, which confirms the renormalizability of the BS 
equation in Eq. fl23|) . 

In fact, we could generalize our approach presented above. One can 
use M in Eq. ( !23l) as the kernel of another BS equation and resum the 
kernel to infinite order once more. Then we would obtain three hierarchies 
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Figure 2: Diagrams included in the second term on the right hand side of Eq. (pO|) . 



of the BS equations. In the same way, these BS equations can be proved 
to renormahzable. Certainly, this procedure can go on to obtain four, five 
hierarchies of the BS equations and so on. This is the hierarchy structure of 
the BS equation. 



3. Propagator and Effective Action 

We have obtained the four-point vertices in the section above, and the 
approximations are made on the level of the four-point functions. Then the 
next task is to find the corresponding two-point function, i.e., the propa- 
gator which constitutes the self-consistent equations, together with the BS 
equations described above. In the following, we will employ the exact renor- 
malization group to obtain the propagator. 

Employing the effective action in Eq. ( ITOj) . taking functional derivative 
with respect to the propagator and considering the stationary condition, we 
have 



5G 

which is in fact the gap equation as follows 



= 0, (31) 



— G'ok ~ ^[Gk], (32) 



9 




where the self-energy S is given by 



6G 



G — G*K 



(33) 



Taking derivative with respect to the flow parameter k on both sides of the 
gap equation in Eq. fj32|) . one arrives at 



The first term is 

The second term can be re-expressed as 



(34) 
(35) 



SG 



dG, 
g=Gk Ok 
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-2i 



SO 



6G5G 



8G~^ 

K 7{ '-^K 

8G~^ 

K 7\ 



G=G^ 



Ok 



1 dG~^ 



Substituting Eqs. ( 1351) ( l36l) into Eq. ( IMI) . one finds 



(36) 



d^G-' = id^R, + -K[G,]Gld,G-\ 



(37) 



Employing the BS equation in Eq. 
arrive at 

d.G-' 



131) . we can solve d^G^^ in Eq. ( 137|) and 



(38) 



+ -M[G^]Gld^R^, 



which is the flow equation for the propagator. Blaizot et al derived this 
equation firstly and found that this fiow equation is completely equivalent 
with the gap equation 
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We should emphasize that in the conventional 
approaches where the approximations are made in the level of the effective 
potential, the gap equation is easily obtained, so the flow equation for the 
propagator seems to be not important. However, in our case the kernel of the 
BS equation is resummed to inflnite order through another BS equation, and 
it is difficult to obtain the gap equation directly. Therefore, we have to resort 
to the fiow equation in Eq. ( 138|) . Then the BS equations ( IT8|) . ( 120|) . (!23|) . and 
the differential equation (138|) constitute a closed system, which can be solved. 

Since the propagator is obtained, the effective action can also be found. 
Differentiating the effective action in Eq. ( ITOl) with respect to the flow pa- 
rameter K, one flnds 



dj:.[G^] 



dT^[G] 



dr^[G] 



G=G^ 



G=Gk 



+ 



dT^[G] 



dG 



dG^ 

G=G^ dn 



Ok 

-]{lT[{d^R^)G^] 



(39) 



where we have used Eq. fl3Tl) in the second line. This is the flow equation for 
the effective action flrst derived in Ref. 
one can get the effective action. 



18| . Integrating the flow equation. 



11 



4. Summary and Outlook 

In this work, we have found the hierarchy structure of the BS equations 
under the formahsm of 2PI effective action theory, and we also employ this 
structure to greatly improve on the approximations for the BS kernel. Re- 
summing the kernel in t and u channels to infinite order is equivalent to 
truncate the effective action to infinite order. Furthermore, our approxi- 
mation approaches do not violate the renormalizability of the theory, which 
is very important for numerical calculations. We also obtain the two-point 
function from the four-point one through fiow evolution equations. There- 
fore, BS equations of different hierarchies and the fiow evolution equation for 
the propagator constitute a closed system, which can be solved completely. 

We should note that there are two hierarchies of BS equations in our 
approach. Therefore, comparing the conventional BS equations, we need 
more computer time in our calculations. But because of remarkable increases 
in computer power made possible by clusters, this problem is not difficult to 
solved. 

Since we have proposed to employ the hierarchy structure of the BS equa- 
tions to improve on the approximations for the BS kernel, it is very interest- 
ing to apply our approaches to detailed problems. One potential interesting 
problem is to compute the shear viscosity of the quark-gluon plasma or other 
thermal fields. As we know the shear viscosity is very sensitive to the proper- 
ties of the four-point vertex, so it is expected that our approach will advance 
the computation of the shear viscosity. Furthermore, our approaches can 
also be applied to strongly correlated electron systems, for example, our ap- 
proaches can be used to study the two-dimensional Hubbard model, which 
may shed new light on our understanding about high superconductors. 
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